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We present a simple discretization by radial basis functions for the Poisson equation 
with Dirichlet boundary condition. A Lagrangian multiplier using piecewise polynomials 
i— i is used to accommodate the boundary condition. This simplifies previous attempts to 

use radial basis functions in the interior domain to approximate the solution and on the 
boundary to approximate the multiplier, which technically requires that the mesh norm 
in the interior domain is significantly smaller than that on the boundary. Numerical 
experiments confirm theoretical results. 
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1 Introduction 

On 

Radial basis functions (RBFs) have been used successfully [11] to solve partial differential 
equations with Neumann or Robin boundary conditions. When Dirichlet conditions are 
considered they must be approximated in an appropriate way. In the case of RBFs this 
causes analytical difficulties since traces of such functions have representations that change 
continuously with the position of their centers. Therefore, corresponding discrete spaces 
(on the boundary) do not have a clear structure that could be used for analysis. Such 
J> difficulties do not appear when treating natural boundary conditions. 

Lagrangian multipliers consisting of RBFs have been studied in [JJ. Due to a weaker 
inverse property of RBFs as compared to that of finite elements, the condition imposed 



C3 



on the mesh norms in the interior domain and on the boundary is too restrictive. In this 
paper we suggest to use RBFs in the domain and finite elements on the boundary. We also 
analyze the influence of scaling of RBFs on the error estimate. Scaling of RBFs can be used 
to avoid too much overlap which is essential for the conditioning of stiffness matrices. 

The idea to couple RBFs with finite elements for the Lagrangian multiplier is proposed 
in [7] where we solve a hypersingular integral equation. In that paper, in order to deal with 
the property that the solution of the integral equation can be extended by zero, we have 
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to consider the Lagrangian multiplier on an extended domain. In the current situation it 
suffices to use multipliers on the boundary of the domain where the problem is set. Moreover, 
not as in [7J where the Lagrangian multiplier does not have any physical meaning, in the 
problem to be considered in this paper it represents the normal derivative of the solution. We 
prove an inf-sup condition which allows an estimate for the approximation of this multiplier; 
we did not succeed in proving this condition for the problem considered in [7J. 

The paper is organized as follows. In Section [2] we introduce the model problem and the 
mixed variational formulation. In Section [3] the finite-dimensional spaces and discretization 
are introduced. Section [4] is a revisit of approximation properties of RBFs where we extend 
previous results so that they can be used in the current study. Section [5] presents the main 
result of the paper, namely an a priori error estimate for the approximation. This analysis 
is carried out after we prove discrete ellipticity of the Dirichlet bilinear form (in the case 
that no mass term is present) and a discrete inf-sup condition of the bilinear form involving 
the Lagrangian multiplier. Section [6] presents numerical experiments that confirm the error 
estimate. 

Throughout the paper, the notation a < b indicates that there exists a constant C > 
independent of discretization or scaling parameters hx, k, r and involved functions (except 
where otherwise noted) such that a < Cb. Similarly we use a > b, and a ~ b means that 



2 Model problem and mixed formulation 

Consider the model problem 

—Au + ku = f in Q, 

(2.1) 

u = g on r, 

where Q C M. d (d = 2, 3) is a polygonal (d = 2) or polyhedral (d = 3) domain with boundary 
r, k > is a constant, and where / G ^(^) an d g G H l l 2 (T) are given functions. When 
k = then / and g must satisfy the usual compatibility condition. 

The solution u of (2.1) will be found in the weak sense by using a standard mixed 
variational formulation. Defining 



a(v, w) 



/ (V v ■ V w + kvw) dx Vu, w G H l (VL), 
Jn 

b( v ,n) = J viids Vv G H l (n), n G H- l/2 (T), 

we can easily see that there exists a constant C > such that 

a{v,w) < C\\v\\ H i in) \\w\\ H i^ \/v,w G H 1 ^), 
b(v,fi) < C\\v\\ Him \\fi\\ H - 1/2(T) Vv G H\n), n G H- l ' 2 {T). 



(2.2) 



(2.3) 



A variational formulation of (2.1) is formulated as: Find (u, A) G i7 1 (fi) x H 1 ^ 2 (T) satis- 
fying 

a(u, v) + b(v, A) = F(v) Vv € -ff 1 (0), 

b(u,n) = G(fj,) V/i G £T 1/2 (r), U '" 
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where 



F(v) = / fvdx and G(/u) = / g^ds. 
Jn Jr 

The following result is well known. 
Proposition 2.1. There exists a unique solution (u, A) G x i7~ 1 / 2 (r) to i/ie probl 



em 



(2.4). Moreover, A = du/dn where n is the outward normal vector on T. 



3 Discretization with RBFs and finite elements 



We first define the finite-dimensional space that approximates u in (2.4). Let $ : M — > M 
be a radial basis function whose Fourier transform $ satisfies 

~ (1 + H 2 )' T , (3.1) 

where r > d/2. We consider the scaled radial basis functions 

$ r (ac) := r" d $(a;/r), r > 0, a; € M d , 

so that 

S P (w) ~ (l + r 2 |a;| 2 )- r , w G K d . (3.2) 
The native space associated with <£ r is defined by 

1 jR d $ r (tt)) J 



where T>'(K d ) is the space of distributions defined in M d . This space is equipped with an 
inner product and a norm defined by 



=/ g( ^ M da, and |Hk := ( [ B^-dJ)^. 
Jm d 3v(w) ^Jm d $r(^) ' 



Under the assumption (3.1 ), the native space A/$ r is isomorphic to the Sobolev space 
with equivalent norm ||v(-)(l + r 2 | - 1 2 ) 9_//2 1| ^ 2 ^d) . 

Given a set of quasi- uniform centers X = {x\, . . . ,xn} C O with mesh norm hx = 
su Pccef2 m ^ n i<j<N \\ x ~ x j\\2, we define 

H x , r = span{$i,...,$iv}, (3.3) 

where $i(a;) = <3? r (:E — ajj) for i = 1, ... , iV. (Note that since the nodes can be near to 
or even on the boundary r of f2, the supports of the scaled radial basis functions are not 



necessarily subsets of f2.) The solution u to (2.4) will be approximated by ux G H X r- 



For the approximation of the Lagrangian multiplier A we use functions (not necessarily 
continuous) which are piecewise polynomials of degree p > defined on a quasi-uniform 
partition Tk of the boundary r of O: 

A fc := {fi : T -> M | /i| T G P p VT G T fc }. (3.4) 
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Here, k is the mesh size of 71, and P p is the space of polynomials of degree at most p. 
Using these discrete spaces, the Galerkin scheme with radial basis functions and La- 



grangian multipliers for the approximate solution of (2.4) is: find (ux,Xk) G Hx, r x 
satisfying 

a(u x ,v) + b(v, A*;) = F(v) Vv £ H x , r , ,„ 

(o.oj 

b(u x ,fj) = G(n) V/x G A fe . 
4 Approximation property of scaled RBFs 

For any integer m > and real r > we denote the norm of the scaled Sobolev space 
H?(Q) by 

, |a|<m 

with multi-index a = {a\, . . . , aj) and \a\ = a% + ■ • • + ad- We note that H®(Q) = H°(Q). 
For s G (0, m) with integer m > we define the scaled interpolation space 

where 9 = s/m. Here we employ the so-called real K- method; see [2J. The interpolation 
norm in H^(Q) can be represented as follows. Let S r ■ Dom(S r r ) C H™(Q) — > H°(tt) be an 
unbounded linear operator defined by 

{S r v,w) H orty = {v,w) Hm ,ty \/v G Dom(SV), Vu> G H™(Q). 

1/2 

It is clear that S r is self-adjoint and positive. Thus there exists A r := S r : H^yVL) — > 
H°(Q) satisfying 

(v,w) H rn( n) = (A r V,A r w} H0{Q) . 

The inner product and norm in H^(Q) can now be represented as 

( v ' w )h^Q) = \ A r V > A r W ) H o, n) alld IMItf^n) = H A r u lliY°(«)- (4-1) 

The Sobolev spaces H*(£l) form a Hilbert scale with the following property. 

Lemma 4.1. Let s% and S2 be non-negative real numbers, and let sq = (s\ + S2)/2. Then 
for any v G H^(U) n H° (Q) there holds 

wev(n)\{o} \\ w \\Hr 2 (n) 

where T>(Q) is the space of all functions which are restrictions on fi of infinitely differentiate 
functions in M. d . 
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Proof. Let m be an integer not less than max{si, S2/- We may assume that maxjsi, 82} > 
0. Then 

where 9i = Sj/m, i = 0, 1, 2. For any v G H^(U) D H^°(Q), there holds 

(^ z >// r s i(n) 



IMIff r s i(fi) = SU P 
zex>(n)\- 

For each 2 G I?(f2) we define if = Af 1 ^ 92 ^. Then by noting (4.1), the relation #0 



(ui -f #2) /2, and the self-adjointness of A^, we obtain 

IkllflJ^n) = l|A^«||flt)(n) = l|A£ 2 Hlno (n) = |Mlflj?(n) 

and 

Therefore, 

\\ v \\H^(n) = SU _P -jj-jj 

wev(n)\{o} \\ w \\Hr 2 (n) 

and the lemma is proved. □ 
For any v£ H^(fl), let Ixv denote its interpolant in the space Hx,r, i-e., ix^ £ #X,r 
satisfies 

Ixv(xj) = v(xj), j = l,...,N. 

For a domain in two dimensions (d = 2), the following approximation property is proved in 
[3 Lemma 5.3]; see also the proof of this lemma. The same arguments apply also to the 
case d = 3. 



Lemma 4.2. Let assumption (3.1) be satisfied. Then for any v G H^{Vt) there holds 



\\v - Ixv\\Hf(a) < C(s,t) \\ v - I x v \\ H -(Q)<C(s,r)( h ^] \\v\\ H r {n) 

for < s < L T J < r < ro wrai/i ro > arbitrary but fixed. 

When v is smoother, the error bound can be extended by using the technique developed 
in [8], [12], and modified in [9] for a sphere. 



Lemma 4.3. Assume that (3.1) holds. Let T be an operator defined by 



Tip(x) := \ $ r (aj - y)ip(y) dy, x G 

Then for any s£R 
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(i) T is an isomorphism from T (R d ) onto H^ +T (M. d ) and satisfies 

ll T ^lli^+-(K<*) - ll^llHr T (R d )' 

(ii) For any tp G H~ T (R d ) and £ G 7V$ r there holds 

(Tip,0$ r = (^,£)i*p(R««) , 
i.e., T is the adjoint of the embedding operator of the native space A/$ r into H®(M. d ). 
Proof. The lemma follows from 

Tty(u) = $ r (u)tp(u)) ~ (1 + r 2 \u\ 2 )- T ^(uj), lj G R d . 

□ 

With the help of the above lemma, we now extend the error bound in Lemma 4.2 when 
v is smoother. 



Lemma 4.4. Let the assumption (3.1) be satisfied. Then for any s,t G M. satisfying < 
s < t < t < 2t, if v £ H*(ri) then the following estimate holds 



(h \ 



V 



for < r < ro with ro > arbitrary but fixed. 

Proof. Consider first the case when t = 2r. Let 

E : H°(£l) -> H?{R d ) 

be an r-uniformly bounded extension operator for any a > 0; cf. [ZJ Lemma 5.1]. For any 
v G H 2r (ft), since IxEv = I x v = EIxv on il, one finds that Ev — IxEv is an extension of 
v — Ixv. Therefore, the property that Ix is an orthogonal projection in A$ T . yields 

\\v - Ixv\\ 2 H r (a) < \\Ev - I x Ev\\ 2 HT(Rd) ~ \\Ev - I x Ev\\% 



< (Ev-I x Ev,Ev)$ 



It follows from Lemma 4.3 that there exists ifi G H®(R d ) such that 

Tip = Ev, \\ip\\ H o {Rd) ~ \\Ev\\ H 2 T{Rd) 

and 

(Ev - IxEv, Ev)$ t = (Ev - I x Ev, tp) H o^ ■ 

Therefore 

\\ v ~ I xv\\ 2 H r {n) < \\Ev - IxEvWjjo^WEvW^r^d) 

< \\v - Ixv\\ H o {n) \\v\\ H 2r {n) . (4.2) 
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By applying Lemma 4.2 with s = and using the above inequality we obtain 

2r 



\v — Ixv\ 



<C(r) 



<C{t) 



hx 
r 

h x 



\v — Ixv\ 



2t 



\v - Ixv\\ H °(a)\\v\\H?T(n)- 



Thus the required estimate is proved for s = and t = 2r. It also follows from (4.2) that 
the required estimate holds for s = r and t = 2r. By using interpolation we deduce the 
estimate for < s < r and t = 2r. Interpolation between 

2t-s 



\v — Ixv\ 



< 



and 



l« - ^x'yllHs(Q) 



< 



X 



I X 



\ v \\H$T(p,) 



\ V \\HZ{Q.) 



□ 



yields the required estimate for 0<s<r<t<2r, proving the lemma. 

To extend Lemma 4.4 to include less smooth functions, i.e. v G H^(Q) for < t < r 
we use Lemma 14.11 



Lemma 4.5. Assume that (3.1) holds. For any v G H^Vt) with < t < 2r i/iere exists 
zx G ifx r satisfying 



v ~ zx\\h«{Q) < C(s,t,r) 



r 



(4.3) 



/or < s < min{i, r} and < r < rp u>ii/i ro > arbitrary but fixed. 



Proof. We only need to prove the lemma for < t < r. Consider first the case when 
r/2 < t < t. Then < 2t - r < t. Let P t : — > H x , r be the projection defined by 



(Ptv,z) 



H*(Q) 



(v,z) m 



fl*(0) 



Vz G %. 



It can be seen that 

\\Ptv - w||fl*(n) < NIb*(«)- 
If s G [0, 2t — t] then r < 2t — s < 2r. Hence, for any w G H^ s it follows from Lemma 

t-s 



(4.4) 
(4.5) 



that 



4.4 



\w - I x w\\ H t An) < 



hx 
r 



(4.6) 



By using Lemma 4.1 and (4.4)-(4.6) we obtain 

\\ p tv - v\\ H s {n) = sup 

wev(tt)\{0} 

= sup 

wev(n)\{o} 

t-s 



(Ptv-v,w) H t {n) 
{P t v -v,w- Ixw) H t 



< 



hx 
r 
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In particular, there holds 

T-t 



\\P t v - v\\ H 2t- T{n) < \\v\\ H t (n) . 



Hence, for s € [2t — T,t] by noting (4.5) and using interpolation we obtain the required 
estimate, and thus prove (4.3) for < s < t and r/2 < t < r. 

By successively considering the case r/4 < t < r/2, then r/8 < t < t/4, etc., and using 
the same argument, we finish the proof of the lemma. □ 

We are now able to derive the approximation property of Hx, r in non-scaled Sobolev 
norms. 



Lemma 4.6. Assume that (3.1) holds. For any v £ iT*(Q) with < t < 2r there exists 
zx 6 Hx,r satisfying 

ht 



\v - z x \\H'(a) < Cis^^Y-^WvWHt^n) 



for < s < min{t, r} and < r < r^ with ro > arbitrary but fixed. 

Proof. First we note that since r G (0, ro], for any function v and any positive integer 
m there hold 

IMI# (C) = IMItfO(n) and ||v||H™(fi)<;iHltf™(fi) < r _m ||t!||^m (Q) . 
By interpolation it follows that for < s < m 

\\Ah°{q)S\\v\\h>>{V) < r~ 8 \\v\\ H a(ny 



The required result is then a consequence of the above inequalities and Lemma 4.5 □ 



5 Error estimate 



In this section we prove our main result, Theorem 5.3, which establishes quasi-optimal con- 
vergence of our mixed method and convergence orders. Of course, approximation properties 
depend on the regularity of solutions. To keep things simple, we assume that, for smooth 
data, we have standard elliptic regularity limited by the smoothness of T. More precisely, 



let 5 be such that the solution u of (2.1 ) for any k > and any sufficiently smooth data /, 
g, satisfies 

S 6(1,2]: u£H 5 (Q). (5.1) 

We restrict our considerations to regularity no more than i7 2 (f2) since we are interested 
in non-smooth problems and to simplify results when approximation spaces use piecewise 
polynomials of higher degree. Of course, in two dimensions f2 being a polygon, there holds 

5 



2 if Q is convex, 

2tt/uj if Q is non-convex, 



with uj G (it, 2tt) being the angle of the largest re-entrant corner of the boundary T in case 
of non-convex 0. The characterization of 5 for a polyhedral domain is a bit more involved 
and not given here. 
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We use standard Babuska-Brezzi theory to prove the main result. In order to do so we 
now prove ellipticity of the bilinear form a(-, •) on 

V X ,r ■= {v G H x ,r) b(v, y) = V/i G A fc } 

when k = and an inf-sup condition for the bilinear form £>(■, •). 

Lemma 5.1. For k sufficiently small there holds 

\ v \m(n) Z \Mm(n) G V x , r - 

Proof. The proof is standard (cf. [5J Lemma 4.5]) and is given for convenience of the 
reader. Let v G Vx,r be given and decomposed as v = vq + D with D = J n v dx so 

that j n vq dx = 0. Here, \Q\ denotes the measure of O. There holds 



(5.2) 



Now let n G -^2(r) \ {0} and its L2(r)-projection H/j, onto be given. By duality, a 
standard approximation result and the trace theorem, we find that there holds 

l&(A M)l = - u o, At)| = \b(v, I 1 ~ n /") ~ K v o, A*) I 

^ ll w lli/i/2(r) fcl/2 ll/ i llL2(r) + INIlL 2 (r)IHlL 2 (r) 

^ (fc 1/2 IMIni(n) + INI|jfi(fi))||Mll£ 2 (r)> 



that is, 



I-DI ^ 



sup 



|6(AM)I 



</t 1/2 l 



^ei 2 (r)\{o} HMllL 2 (r) 



Ih 1 ^) + ll^ollifi(n)- 



(5.3) 



Now, using (5.2), (5.3) and Poincare-Friedrichs' inequality, we find that 

|2 



,,2 ^ II ||2 i ; || ||2 , 

U\\uun\ 5, WVoWui,^ + k\\v\\ H l( n ) < |r 



Selecting small enough finishes the proof. 



if!(n) 



+ fcl u| 



ff!(fi)- 



□ 



Lemma 5.2. Assume that (3.1) and i/ie elliptic regularity (5.1) /io/d. Suppose that the 
parameters hx, k and r are chosen such that 



K t (hx, k,r) 



is sufficiently small 



(5.4) 



with t = 5 when 5 < 2, or with some t G [1,2) when 5 = 2. Then there exists a positive 
constant a, independent of h x , k and r, except for their relation via Kt, such that there 
holds 

b(v x ,Hk) 

a ll/ i fcll,H--i/2(r) 



sup 

»x6ffx, r \{o) ll^||^i(n) 



V> fc G A fe . 
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Proof. For any ^ £ A^, consider the problem 

—Aw + w = in U, 

dw _, (5-5) 

— = [i k on r. 

on 

Since u k G H S (F) for any s < 1/2, there exists a unique variational solution w G if (fi) of 
(5.5) with regularity estimate (limited by assumption (5.1) depending on the geometry) 

||w||fftnta{«.3/2+*}(fi) ^ |lMfe||ffmta{<5-3/2,a}(r) Vs G [-1/2,1/2); (5.6) 

see e.g. [6]. Moreover, it is shown in [1, Theorem 2.7] that 



b(w,u k ) ~ ||^fc||^-i/2 (r) - 



(5.7) 



On the other hand, since 1 < 5 < 2 < 2r we can invoke Lemma 4.6 to obtain, for some 

ZX G Hx,r, 

with 



\w - zxWhho) < c{t)—r\\M\H*(n) 



t = 5 if<5<2, 
Vt€[l,2) if 5 = 2. 



By using (5.6), the inverse property 



and the assumption (5.4), we deduce that 



t-i 



ti 

\ w ~ z x\\m(Q) < r ttf-i ll^fcllff-i/^r) = K t{hx-,k,r)\\u k \\ H -i/2 { Y^ 



(5.8) 



and (5.6) with s = —1/2 give 



with t as given (if 5 = 2) or chosen (if 8 < 2) for ( |5.4| ). This inequality, assumption ( |5.4| 

(5.9) 



On the other hand (2.3), (5.7) and (5.8) yield for Kf(h x , k,r) small enough 
b(z x ,u k ) = b(w,u k ) + b(z x - w,n k ) 

~ 1 1 ^fc I la-- l/a ( r ) - ll^fc||ff-i/2(r)||^ - ^x||ffi(n) 
> (1 - K t (h x ,k,r)) \\u k \\ 2 H _ 1/2{r) > ||Mfc||ff-i/2( r ), 



i.e., 



\fi>k\\H-V 2 (r) 



< 



b(z x ,u k ) 
l^fc||ff-i/2(r) 



This together with (5.9) yields 



b(v x ,u k ) 



b(z x ,u k ) 

|Mfe||ff-V2(r) ^ I, — n < SU P „ „ 

FX || ff i(O) vx&Hx r\{0} H^Xllffi(fJ) 



□ 



proving the lemma. 

The following theorem is our main result. It proves the quasi-optimal convergence of 
our mixed RBF approximation of the Dirichlet problem with finite element Lagrangian 
multiplier. 
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Theorem 5.3. Let us assume that (3.1) and (5.4) hold and consider radius parameters 
r > which are bounded. In the case that k = in (2.1), we additionally assume that k 
is small enough. Then there exists a unique solution (ux,\k) £ Hx,r x Afc to £/ie problem 
(3.5). Moreover, let f and g be sufficiently smooth so that (u, A) is £/ie solution to (2.4) 
with u G # a (0) and A G If 5-3 / 2 (r) tf G (1,2], cf. Kty. TTien 



^xllifi(n) + ||A - Afe|| ff -i 



/2(D £ 



inf 



^-1 _ 

~ ~^5~W u \\H s (n) + & ll^lljy*-3/2(r)> 



where the implicitly appearing constants depend on the constant C in (2.3), a in Lemma 5.2, 
and the ellipticity constant of bilinear form a(-, •). 

Proof. For k > 0, the first bound of quasi-optimal convergence follows from standard 
Babuska-Brezzi theory (cf. [3] Corollary 12.5.18]) by making use of the continuity and 
ellipticity of the bilinear form a(-, •), and continuity and continuous as well as discrete 



inf-sup condition of the bilinear form b(-, •), cf. Lemma 5.2 
In the case that k 



we use Lemma 5.1 to obtain V3c |r -ellipticity of a(v,v) 



\ Vv\ 2 dx. Then the result follows the same way. 

To show the second estimate for k > we use Lemma 4.6 and the approximation 
property of □ 
Depending on the regularity of the solution, the error estimate above can be simplified. 



Corollary 5.4. Let the assumptions of Theorem 5.3 be satisfied. If 5 > 3/2 then 

,5-1 

\ u \\H s (n)- 



\ u ~ u x\\m(n) + P - ^fclln-vajT) ~ ( 



X 



+ k 



Proof. In the case u G H S (Q) with S > 3/2, the normal derivative A = du/dn can be 
defined in the standard way (normal component trace of weak gradient) so that 



The assertion then is direct consequence of Theorem 5.3 



□ 



6 Numerical results 



We consider the model problem (2.1) with SI = (0, 1) x (0, 1), k = and /, g such that 
u(x\,X2) = x\ + £c|, he. u G H 5 {Q) with 5 = 2. The nodes of X are distributed uniformly 
on including nodes on the boundary. 

We use scaled radial basis functions with the radial basis functions defined in [10| . We 
consider the two cases r = 1.5 and r = 2.5 which correspond to C° and C 2 -functions, 
respectively. They are rotations of univariate polynomials of degrees 2 and 5, respectively. 
We have implemented the method by numerical integration with an overkill of number of 
integration nodes. 

With 5 = 2 assumption (5.4) requires that the ratio fc i_-f r 2- e be small enough for 



some e > 0. We simply choose k ~ hx, r / r (more precisely an integer approximation to 
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1/k smaller than or equal to 1 since the length of the sides of is 1). In this way, for 
fixed r, Kt is fixed and ( |5.4[ ) is not guaranteed. However, our numerical results do not 
show stability problems (that might be caused by a violation of the inf-sup condition) in 
the range of unknowns under consideration. 

For fixed r, the error estimate derived in Corollary |5.4| gives an upper bound 

\\ u - u x\\m(n) + 11^ _ ^fcll,H--i/2(r) <hx + k. 

In the graphs below we plot the individual errors ||u — uxWh 1 ^) an d ||A — Afc||^ 2 (p) on a 
double logarithmic scale. For the latter error, which is measured in the L<i rather than 
i7 _1 / 2 -norm, we expect a reduced convergence like /c 1 / 2 . Both expected error terms, hx 
(labeled as h) and A; 1 / 2 , are also given in the plots (multiplied by 10 to shift them closer to 
the corresponding error curves). Figures [T] and [2] show the results for r = 0.2 with r = 1.5 
and r = 2.5, respectively. Figures [3] and [4] show the corresponding results for reduced 
radius r = 0.1. In all the cases there is some pre- asymptotic range and the errors behave 
as expected for larger number of unknowns. 




Figure 1: Errors for r = 0.2 and r = 1.5. 
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Figure 3: Errors for r = 0.1 and r = 1.5. 




Figure 4: Errors for r = 0.1 and r = 2.5. 
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